Quantum Convolutional Coding with Shared 
Entanglement: General Structure 



Mark M. Wilde and Todd A. Brun 



OO 

O 

o 

(N 



^1" 
(N 



43 ■ 



> 

m 

o 

en 

o 
oo 
o 



X 



Abstract — We present a general theory of entanglement- 
assisted quantum convolutional coding. The codes have a convo- 
lutional or memory structure, they assume that the sender and 
receiver share noiseless entanglement prior to quantum commu- 
nication, and they are not restricted to possess the Calderbank- 
Shor-Steane structure as in previous work. We provide two 
significant advances for quantum convolutional coding theory. We 
first show how to "expand" a given set of quantum convolutional 
generators. This expansion step acts as a preprocessor for a 
polynomial symplectic Gram-Schmidt orthogonalization proce- 
dure that simplifies the commutation relations of the expanded 
generators to be the same as those of entangled Bell states 
(ebits) and ancilla qubits. The above two steps produce a set of 
generators with equivalent error-correcting properties to those 
of the original generators. We then demonstrate how to perform 
online encoding and decoding for a stream of information qubits, 
halves of ebits, and ancilla qubits. The upshot of our theory is that 
the quantum code designer can engineer quantum convolutional 
codes with desirable error-correcting properties without having 
to worry about the commutation relations of these generators. 

Index Terms — quantum convolutional codes, entanglement- 
assisted quantum convolutional codes, quantum information 
theory, entanglement-assisted quantum codes 

I. Introduction 

The theory of quantum convolutional coding has recently 
emerged as a powerful means for protecting a stream of quan- 
tum information from the negative effects of a noisy quantum 
communication channel [Q, El, 0, G), 0, 0, 0, 0, 0, 
iflOl . ifTTI . II 1 21 . We believe that quantum convolutional coding 
theory is one step along the way to finding quantum error- 
correcting codes that approach the capacity of a noisy quantum 
communication channel for sending quantum information 1 13 1, 
lfT4l . Ifl5l |[T6l . H3, ED, GE. Poulin et al. have recently 
incorporated some of this well-developed theory of quantum 
convolutional coding into a theory of quantum serial-turbo 
coding 1 20] with the goal of designing quantum codes that 
come close to achieving capacity. 

Quantum convolutional codes have numerous benefits. The 
periodic structure of their encoding and decoding circuits 
ensures a low complexity for encoding and decoding while 
also providing higher performance than a block code with 
equivalent encoding complexity [7|. The encoding and de- 
coding circuits have the property that the sender Alice and 
the receiver Bob can respectively send and receive qubits in 
an "online" fashion. Alice can encode an arbitrary number of 
information qubits without worrying beforehand how many 
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she may want to send over the quantum communication 
channel. 

A seminal paper on quantum error correction showed how 
to produce quantum error-correcting codes from classical 
linear binary or quaternary block codes 11211 . Importing binary 
codes to construct quantum codes is now known as the 
Calderbank-Shor-Steane (CSS) construction 12211 . Importing 
classical codes is useful because quantum code designers 
can utilize classical block codes with high performance to 
construct quantum codes with high performance. The only 
problem with this technique is that the imported classical block 
codes have to satisfy a restrictive dual-containing condition so 
that the operators in the resulting quantum code commute with 
one another. 

Brun, Devetak, and Hsieh alleviated the dual-containing 
constraint by assuming that the sender and receiver share 
noiseless entanglement in the entanglement-assisted stabi- 
lizer formalism 11231 . l24l . The entanglement-assisted stabi- 
lizer formalism allows quantum code designers to produce 
a quantum code from an arbitrary classical binary or qua- 
ternary block code by incorporating shared entanglement. In 
a different context, Forney et al. showed how to produce 
quantum convolutional codes from classical convolutional 
codes by extending the ideas from the CSS construction 
to the convolutional setting (7); but again, these imported 
classical convolutional codes have to satisfy the restrictive 
dual-containing constraint in order to form valid quantum 
codes (the dual-containing constraint is actually quite a bit 
more restrictive in the convolutional case). The current authors 
followed these two developments by detailing entanglement- 
assisted quantum convolutional codes [11] that have the CSS 
structure. We demonstrated in Ref. ifTll how to produce an 
entanglement-assisted quantum convolutional code from two 
arbitrary classical binary convolutional codes — these imported 
codes do not need to satisfy a dual-containing constraint. The 
benefit of all these entanglement-assisted techniques is that 
we can use high-performance classical codes for correction of 
quantum errors. 

In this paper, we show how to encode and decode an 
entanglement-assisted quantum convolutional code that does 
not necessarily have the CSS structure. The methods in this 
work represent a significant extension of our previous work 
on CSS entanglement-assisted quantum convolutional codes 
Him . In particular, we develop a set of techniques that make 
the commutation relations for the generators of a general code 
the same as those of entangled qubits (ebits) and ancilla qubits. 
This procedure first "expands" the check matrix for a quantum 
convolutional code and applies an extended version of the 
symplectic Gram-Schmidt orthogonalization procedure 11241 



that incorporates the binary polynomials representing quantum 
convolutional codes. We then show how to encode a stream of 
information qubits, ancilla qubits, and ebits to have the error- 
correcting properties of the desired code. We follow by detail- 
ing the decoding circuit that Bob employs at the receiving end 
of the channel. The algorithms for encoding and decoding use 
similar techniques to those outlined in Refs. J3], Il25l . ATI . The 
encoding circuits incorporate both finite-depth and infinite- 
depth ifTTl operations and the decoding circuits incorporate 
finite-depth operations only. One benefit of the techniques 
developed in this paper is that the quantum code designer 
can produce an entanglement-assisted quantum convolutional 
code from a classical quaternary convolutional code. More 
generally, quantum convolutional code designers now have the 
freedom to design quantum convolutional codes with desirable 
error-correcting properties without having to search for codes 
that satisfy a commutativity constraint. 

We structure our work as follows. Section HI] reviews the 
general definition of a quantum convolutional code and the ma- 
trix of binary polynomials that represents it. It also discusses 
the shifted symplectic product that captures the commutation 
relations of Pauli sequences. Section ITTT1 reviews the theory of 
entanglement-assisted quantum block coding. In Section [TV] 
we present an example that demonstrates how to expand the 
check matrix of a set of quantum convolutional generators and 
then show how to expand an arbitrary check matrix. We show 
in Section |V] how to compute the symplectic Gram-Schmidt 
orthogonalization procedure for the example in Section [IV] 
and then generalize this procedure to work for an arbitrary 
set of quantum convolutional generators. The Gram-Schmidt 
orthogonalization technique reduces the quantum check matrix 
to have the same commutation relations as a set of ebits and 
ancilla qubits. This technique is essential for determining an 
encoding circuit that encodes a stream of information qubits, 
ancilla qubits, and ebits. Section [VT] gives the algorithms for 
computing the encoding and decoding circuits for an arbi- 
trary entanglement-assisted quantum convolutional code. We 
present a detailed example in Section [VTTI that illustrates all of 
the procedures in this article. We finish with a discussion of the 
practical issues involved in using these entanglement-assisted 
quantum convolutional codes and present some concluding 
remarks. 

II. Review of Quantum Convolutional Codes 

We first review the theory of quantum convolutional coding. 
The reader can find more detailed reviews in Refs. 0, iflOll . 
IfTTl . A quantum convolutional code admits a representation 
in terms of Pauli sequences — infinite tensor products of Pauli 
matrices. We demonstrate how to map these Pauli sequences 
to a finite-dimensional binary polynomial matrix and call this 
mapping the Pauli-to-binary (P2B) isomorphism. It is more 
straightforward to work with binary polynomial matrices rather 
than Pauli sequences because we can employ linear-algebraic 
techniques. 

We state some basic definitions and notation before pro- 
ceeding with the definition of a quantum convolutional code. 
A Pauli sequence A is a countably infinite tensor product of 



Pauli matrices: 



)A, 



i=0 

where each A4 is a matrix in the Pauli group II = {/, X, Y, Z}. 
Let II denote the set of all Pauli sequences. A finite-weight 
Pauli sequence is one in which the operators X, Y, and Z 
act on a finite number of qubits and the identity operator acts 
on the countably infinite remaining qubits. An infinite- weight 
sequence is one that is not finite-weight. 

Definition 1: A rate-k/n quantum convolutional code ad- 
mits a representation by a basic set Qo of n — k generators 
and all of their n-qubit shifts. For a quantum convolutional 
code, the generators commute with each other and all of the 
n-qubit shifts of themselves and the other generators. This 
commutativity requirement is necessary for the same reason 
that standard stabilizer codes require it l22l . The generators 
in the basic set Qo are as follows: 

Qo = JGi G n z+ : 1 < i < n - k\ . 

A frame of the code consists of n qubits so that the frame 
size is n. The definition of a quantum convolutional code as 
n-qubit shifts of the basic set Qo is what gives the code its 
periodic structure. 

We do not discuss the detailed operation of a quantum 
convolutional code in this article. Please consult Refs. Q, 
irTol . IfTTl for a discussion of a quantum convolutional code's 
operation and see Figure 3 of Ref. iflOl for a depiction of its 
operation. The current article addresses the issues of quantum 
convolutional code design. 

We illustrate the P2B isomorphism by means of an example. 
Consider the following generator for a quantum convolutional 
code: 

• • • | III | XXX | XZY | III | • • • 

We obtained the above generator by modifying generators for 
a code from Ref. [7|. The vertical bars indicate that the frame 
size is three qubits and every three-qubit shift of the above 
generator yields another stabilizer generator for the code. The 
above generator forms a valid quantum convolutional code 
because it commutes with all of its three-qubit shifts. For now, 
we are not concerned with the error-correcting properties of 
the above generator — we are merely interested in illustrating 
the principle of the P2B isomorphism. 

The above generator has the following representation as a 
binary polynomial matrix: 

[ D D | 1 + D 1 1 + D ]. 

The matrix on the right is the "X" matrix and the matrix on 
the left is the "Z" matrix. The vertical bar in the above matrix 
serves to distinguish between the "X" and "Z" matrices. Each 
"X" operator in the first frame of the code gives a "1" entry 
in the "X" matrix. There are no "1" entries in the "Z" matrix 
because there are no "Z" operators in the first frame of the 
code. Each "X" and "Z" operator in the second frame of the 
code gives a "D" entry in the respective "X" and "Z" matrices. 
The "Y" entry in the second frame contributes a "I?" entry to 
both the "X" and "Z" matrix. 



In general, we represent a rate-fc/n quantum convolutional 
code with an (n — k) x 2n-dimensional quantum check matrix 
H (D) whose entries are binary polynomials where 

H(D) = [ Z(D) | X(D) ]. 

The shifted symplectic product captures the commu- 
tations relations of any two generators ifTUl . ifTTI . Sup- 
pose that hi(D) = [ zi(D) \ xi(D) ] and h 2 (D) = 
[ z 2 (D) I x 2 (D) ] are the first and second respective 
rows of check matrix H (D). The shifted symplectic product 
(hi © h 2 ) (D) is a polynomial equal to the following expres- 



(h x h 2 ) (D) = z x (D) ■ x 2 (D- 1 ) + xi (D) ■ z 2 (D- 1 ) , 

where "•" represents the standard inner product between two 
vectors and addition is binary. Please note that we have slightly 
changed the convention for the shifted symplectic product 
from that used in Ref's. iflOl . IfTTI in order to agree with the 
convention in Ref's. 0, 0. The shifted symplectic products 
(hi hi) (D), (h 2 h 2 ) (D), and (hi h 2 ) (D) are equal to 
zero for our above example because these generators are part 
of a valid quantum convolutional code and therefore commute 
with each other and with all shifts of themselves and each 
other. For a general set of generators, this condition does not 
necessarily have to hold. We spend much of our effort in 
this article developing techniques to modify the commutation 
relations of a general set of generators. 

The following matrix il (D) captures the commutation 
relations of the generators in H (D): 

n(D) = Z (D) X T (D- 1 ) + X (D) Z T (D- 1 ) . 

The reader can verify that the matrix elements [J7 (£>)]■ of 
il (D) are as follows: 

P(D)] ij = (h i Qh j )(D), 

where denotes the shifted symplectic product. The matrix 
n (D) obeys the symmetry: Q (D) = fl T (-D -1 ). We call the 
matrix f2 (D) the shifted symplectic product matrix because 
it encodes all of the shifted symplectic products or the com- 
mutation relations of the code. This matrix is equal to the 
null matrix for a valid quantum convolutional code because 
all generators commute with themselves and with all n-qubit 
shifts of themselves and each other. For a general set of 
generators, il (D) is not necessarily equal to the null matrix. 
We comment briefly on row operations. We can add one row 
of the check matrix H (D) to another row in it without chang- 
ing the error-correcting properties of the code. This property 
holds because the code is additive. We can even premultiply 
check matrix H (D) by an arbitrary matrix R (D) with rational 
polynomial entries without changing the error-correcting prop- 
erties of the code. This result follows from the classical theory 
of convolutional coding 1261 , Let H ' (D) denote the resulting 
check matrix where H ' (D) = R(D)H(D). The resulting 
effect on the shifted symplectic product matrix Q, (D) is to 
change it to another shifted symplectic product matrix Q! (D) 
related to SI (D) by 



Row operations do not change the commutation relations 
of a valid quantum convolutional code because its shifted 
symplectic product matrix is equal to the null matrix. But 
row operations do change the commutation relations of a set 
of generators whose corresponding shifted symplectic product 
matrix is not equal to the null matrix. This ability to change 
the commutation relations through row operations is crucial 
for constructing entanglement-assisted quantum convolutional 
codes from an arbitrary set of generators. We use entanglement 
to resolve any anticommutativity in the generators. 

III. Review of Entanglement-Assisted Quantum 
Block Codes 

We briefly review the theory of entanglement-assisted quan- 
tum coding for block quantum codes. The entanglement- 
assisted stabilizer formalism subsumed the stabilizer formal- 
ism ll27l . II2TI by assuming that the sender and receiver in 
a quantum communication protocol can share entanglement 
[24 1, |23l . This technique of incorporating shared entan- 
glement allows the sender and receiver to import arbitrary 
classical codes for use in quantum error correction. 



Consider the ebit or Bell state |$+) 
sender Alice A and the receiver Bob B 



AB 



shared between the 



|*+) 



AB 



|0)>) B 



ll>' 



V2 



Two operators that stabilize this state are X A X B and Z A Z B . 
These two operators commute, 

[X A X B ,Z A Z B ] =0, 

but consider that the local operators operating only on either 
party anticommute, 

{X A ,Z A } = {X B ,Z B }=0. 

The above commutation relations hint at a way that we can 
resolve anticommutativity in a set of generators. Suppose that 
we have two generators in a code that anticommute. We can re- 
solve the anticommutativity in the two generators by adding an 
extra X 01 Z operator to the generators on the receiver's side, 
at the cost of an ebit of entanglement. We can find an encoding 
circuit starting from our unencoded information qubits and one 
ebit to encode the generators for the code. Detailed algorithms 
exist to find an encoding circuit — please see Refs. [24|, |23|, 
1 10], [28 1 . These algorithms can also handle a set of generators 
that have more complicated commutation relations. It turns out 
that these algorithms use the optimal number of ebits for a 
given set of generators [291 . 

Let us write the local operators X A and Z A as respective 
symplectic vectors gA and h^- 

9 A = [ I 1 ] , 
h A = [ 1 J ] . 

The symplectic product matrix il of these two symplectic 
vectors is as follows: 



Q' (D) =R(D)fl (D) R T (d- 1 ) 



(1) 



n = 



1 

1 



(2) 



The above symplectic product matrix is so special for the 
purposes of this article that we give it the name J = fi. It 
means that two generators have commutation relations that 
are equivalent to half of an ebit that the sender possesses. 

A general set of generators for a quantum block code 
can have complicated commutation relations. There exists 
a symplectic Gram-Schmidt orthogonalization procedure that 
simplifies the commutation relations tf24), ill], JlO), J28). 
Specifically, the algorithm performs row operations that do 
not affect the code's error-correcting properties and thus gives 
a set of generators that form an equivalent code. For a given 
set of n — k generators, the commutation relations for the 
generators resulting from the algorithm have a special form. 
Their symplectic product matrix il has the standard form: 



^ = 0^®0[O], 



(3) 



£=1 



3=1 



where the large and small correspond to the direct sum 
operation and [0] is the one-element null matrix. The standard 
form above implies that the commutation relations of the n—k 
generators are equivalent to those of c halves of ebits and a 
ancilla qubits where a = n — k — 2c, It is straightforward to 
find an encoding circuit that encodes k + c information qubits, 
a ancilla qubits, and c halves of ebits once we have reduced 
the commutation relations of the generators to have the above 
standard form. 

The focus of the next two sections is to elucidate techniques 
for reducing a set of quantum convolutional generators to 
have the standard commutation relations given above. The 
difference between the current techniques and the former 
techniques l24l . Il23l . £l0|, |28| is that the current techniques 
operate on generators that have a convolutional form rather 
than a block form. The next section shows how to expand 
a set of quantum convolutional generators to simplify their 
commutation relations. The section following the next outlines 
a symplectic Gram-Schmidt orthogonalization algorithm that 
uses binary polynomial operations to reduce the commutation 
relations of the expanded generators to have the standard form. 

IV. The Expansion of Quantum Convolutional 
Generators 

We begin this section by demonstrating how to expand 
a particular generator that we eventually incorporate in an 
entanglement-assisted quantum convolutional code. We later 
generalize this example and the expansion technique to an 
arbitrary set of generators. This technique is important for 
determining how to utilize entanglement in the form of ebits 
in a quantum convolutional code. 

A. Example of the Expansion 

Let us first suppose that we have one convolutional gener- 
ator: 



I X 



I 



(4) 



This generator does not yet represent a valid quantum convo- 
lutional code because it anticommutes with a shift of itself by 
one qubit to the left or to the right. We are not concerned with 



the error-correcting properties of this generator but merely 
want to illustrate the technique of expanding it. 

Let us for now consider a block version of this code that 
operates on six physical qubits with six total generators. The 
generators are as follows: 



X 

I 
I 
I 
I 
I 



z 


I 


I 


I 


X 


z 


I 


I 


I 


X 


z 


I 


I 


I 


X 


z 


I 


I 


I 


X 


I 


I 


I 


I 



I 
I 
I 
I 
z 

X 



(5) 



We still use the vertical bars in the above block code to denote 
that the frame size of the code is one. Observe that we can 
view the frame size of the code as two without changing any 
of the error-correcting properties of the block code: 



X 

I 
I 
I 
I 
I 



Z 
X 

I 
I 
I 
I 



I 
z 

X 

I 
I 
I 



I 
I 

z 

X 

I 
I 



I 
I 
I 
z 

X 

I 



I 
I 
I 
I 
z 

X 



(6) 



The frame is merely a way to organize our qubits so that 
we send one frame at a time over the channel after the 
online encoding circuit has finished processing them. We can 
extend the above block code with frame size two to have a 
convolutional structure with the following two convolutional 
generators: 



I I 


X 


z 


I I 


I I 


I I 


I 


X 


Z I 


I I 



The above two generators with frame size two have equivalent 
error-correcting properties to the original generator in © by 
the arguments above. We say that we have expanded the 
original generator by a factor of two or that the above code 
is a two-expanded version of the original generator. We can 
also extend the original generator in (0]i to have a frame size 
of three and we require three convolutional generators so that 
they have equivalent error-correcting properties to the original 
generator: 

III X Z I III III 
III I X Z III III 
III I I X Z I I I I I 

The representation of the original generator in © as a 
quantum check matrix in the polynomial formalism is as 
follows: 

g (D) = [ D | 1 ] . (7) 

The two-expanded check matrix has the following polynomial 
representation: 



G 2 (D) = 



1 
D 



1 
1 



(8) 



and the three-expanded check matrix has the following poly- 
nomial representation: 



G 3 (D) 






1 





1 


D 






1 
1 
1 



An alternative method for obtaining the polynomial repre- 
sentation of the two-expanded check matrix consists of two 
steps. We first multiply g (D) as follows: 



GUD) 



1 
D 



[9(D)] 



1 D- 1 






1 D- 1 



We then "plug in" the fractional delay operator D 1 / 2 and apply 
the flooring operation |_ - J that nulls the coefficients of any 
fractional power of D: 



G 2 (D) = 



G< 2 (D^)\ 



A similar technique applies to find the check matrix of the 
three-expanded matrix. We first multiply g (D) as follows: 



G 3 (D) 



1 
D 

D 2 



[9(D)] 



1 5 F 

1 J- J- 
1 D D 2 







We then "plug in" the fractional delay operator I? 1 / 3 and apply 
the flooring operation |_-J : 



G 3 (D) 



G' 3 (D 1 / 3 



We discuss the general method for expanding an arbitrary 
check matrix in the next subsection. 



B. General Technique for Expansion 

We generalize the above example to determine how to 
expand an arbitrary quantum check matrix by a factor of I. 
Suppose that we have an n — k x 2n quantum check matrix 
H (D) where 

H(D)=[ Z(D) | X(D) ]. 

Let D denote a diagonal matrix whose diagonal entries are 
the delay operator D, We take the convention that D° is the 
identity matrix and D m is the matrix D multiplied m times 
so that its diagonal entries are D m . Let Ri (D) and Ci (D) 
denote the following matrices: 



Rl (D) = 


D° D 1 


D< 


- 1 ] i = 




Ci (D) = 


D° 

••• 





• 
D° • 




• D-C- 1 



where the diagonal D matrices in Ri (D) and Ci (D) have 
respective dimensions n — kxn — k and nxn. We premultiply 
and postmultiply the matrix H (D) by respective matrices 
Ri (D) and Ci (D), evaluate the resulting matrix at a fractional 
power \/l of the delay operator D, and perform the flooring 
operation |_-J to null the coefficients of any fractional power 
of D. The ^-expanded check matrix Hi (D) is as follows: 



,1/1 



)V 



,i/i 



Hi{D)= Ri[D^ l )H[D^ l )Ci[D 



The ^-expanded quantum check matrix Hi (D) has equivalent 
error-correcting properties to the original check matrix. 



V. Polynomial Symplectic Gram-Schmidt 
Procedure 

In general, a given set of generators may have complicated 
commutation relations. We have to simplify the commutation 
relations so that we can encode information qubits with the 
help of ancilla qubits and halves of ebits shared with the 
receiver. In this section, we begin with an arbitrary set of 
convolutional generators. We show how to determine a set 
of generators with equivalent error-correcting properties and 
commutation relations that are the same as those of halves 
of ebits and ancilla qubits. We first show the technique for 
an example by illustrating it for Pauli sequences, for the 
quantum check matrix, and with the shifted symplectic product 
matrix. We then state a polynomial symplectic Gram-Schmidt 
orthogonalization algorithm that performs this action for an 
arbitrary set of quantum convolutional generators. 

A. Example of the Polynomial Symplectic Gram-Schmidt Or- 
thogonalization Procedure 

1) Pauli Picture: Let us consider again our example from 
the previous section. Specifically, consider the respective ex- 
pressions in (0]l and (O for the convolutional generator and the 
block code. Recall that we can multiply the generators in a 
block code without changing the error-correcting properties of 
the code [22]. Therefore, we can multiply the sixth generator 
in (O to the fourth. We can also multiply the sixth and the 
fourth to the second to yield the following equivalent code: 

X 
I 
I 
I 
I 
I 



z 


I 


I 


X 


z 


X 


I 


X 


z 


I 


I 


X 


I 


I 


I 


I 


I 


I 



I 


I 


z 


X 


I 


I 


z 


X 


X 


z 


I 


X 



(9) 



We have manipulated the two-expanded matrix rather than 
the code in © because the commutation relations of the 
above code are equivalent to the commutation relations of the 
following operators: 



z 


X 


I 


I 


I 


I 



I 


I 


I 


I 


I 


z 


I 


X 


I 


I 


I 


I 



I 


I 


I 


I 


z 


X 



We can use three ebits to encode the set of generators in 
(O because they have the same commutation relations as the 
above operators and the above operators correspond to halves 
of three ebits. We resolve the anticommutation relations by 
using the following entanglement-assisted code: 



Z X 
X I 



z 

X 

I 
I 
I 
I 



I 
I 
z 

X 

I 
I 



I I 

Z X 

X z 



X 

I 
I 



I 
I 
I 
I 

Z X 
X I 



I I 

Z X 

I I 

Z X 

z 

X 



The convention above is that the first qubit of each frame 
belongs to Bob and corresponds to half of an ebit. The second 
two qubits of each frame belong to Alice. The overall code 
forms a commuting stabilizer so that it corresponds to a valid 
quantum code. Bob could measure the above operators to 
diagnose errors or he could measure the following operators 
that are equivalent by row operations: 



Z X 



X 
I 
I 
I 
I 



z 

X 

I 
I 
I 
I 



I 

X 

z 

X 

I 
I 



I 
z 

X 

I 
I 
I 



I 
I 

z 

X 

I 
I 



I 
I 
I 

X 

z 

X 



I 
I 
I 
z 

X 

I 



I 
I 
I 
I 

z 

X 



One can check that the operators corresponding to the sec- 
ond two qubits of each frame are equivalent to the desired 
generators in ((6]). 

2) Polynomial Picture: Let us consider the convolutional 
generator in ®. We now use the polynomial formalism 
because it is easier to perform the manipulations for general 
codes in this picture rather than in the Pauli picture. 

We extend the row operations from the above block code 
to the polynomial picture. Each row operation multiplied 
the even-numbered generators by themselves shifted by two 
qubits. Extending this operation to act on an infinite Pauli 
sequence corresponds to multiplying the second generator in 
dHJ by the rational polynomial 1/ (1 + D). Consider the two- 
expanded check matrix with the operation described above 
applied to the second generator: 



" 9i (D) ' 




1 


1 





_ 92 (D) 




-2- 





1 

l + D 



The shifted symplectic products (g\ g\) (D) = 0, 
(gi&g 2 ){D) = 1, and (g 2 Qg 2 )(D) = capture the 
commutation relations of the resulting code for all frames. We 
can therefore resolve this anticommutativity by prepending a 
column that corresponds to Bob possessing an ebit: 





D 

l+D 





1 

l+D 



Bob possesses the qubit corresponding to column one of both 
the "Z" and "X" matrix and Alice possesses the two qubits 
corresponding to the second and third columns. The code 
above has equivalent error-correcting properties to those of 
the desired generators. 

The generators in the above polynomial setting correspond 
to Pauli sequences with infinite weight. Infinite-weight gener- 
ators are undesirable because Bob cannot measure an infinite 
number of qubits. There is a simple solution to this problem 
and it is similar to what we did at the end of the previous 
subsection. We multiply the second row of the above check 
matrix by 1 + D to obtain the following check matrix that has 
equivalent error-correcting properties to the above one: 



1 1 
D 



1 
l + D 1 



The generators in the above check matrix correspond to 
Pauli sequences with finite weight by the P2B isomorphism. 



Specifically, the above two generators are equivalent to the 
following two Pauli sequences and all of their three-qubit 
shifts: 



Z 
X 



X 

I 



z 

X 



I 

X 



I 

z 



Bob can measure these operators because they have finite 
weight. 

3) Shifted Symplectic Product Matrix Picture: We can more 
easily determine the row operations that simplify the com- 
mutation relations by looking only at the shifted symplectic 
product matrix. We would like to perform row operations on 
the check matrix to reduce its corresponding shifted symplectic 
product matrix to the standard form in (0, so that it has 
commutation relations equivalent to those of halves of ebits 
and ancilla qubits. 

The shifted symplectic product matrix corresponding to the 
check matrix in (0 is a one-element matrix fl g (D) where 

fi s (D) = [D- 1 + D] . 

It is clear that fl g (D) is not reducible by row operations to 

the 2x2 matrix J in (O because il g (D) is a one-element 

matrix. It is also not reducible to the one-element null matrix 

[0] because there is no valid row operation that can zero out 

the term D~ Y + D. 

We therefore consider the two-expanded check matrix in ([8]) 

to determine if we can reduce it with row operations to the 

standard form in (|3). The shifted symplectic product matrix 

riG 2 (D) for the two-expanded check matrix G2 (D) is as 

follows: 

1 + D- 

1 + D 



"g 2 (D) 



R{D) 



We can represent the row operation of multiplying the second 
generator by 1/ (1 + D) as a matrix R (D) where 

1 

l/(l + D) 

The effect on the matrix Q,g 2 (D) is to change it to 

R{D)n G2 {D)R T (D- 1 )= J, 

as described in (|TJ- The above matrix J has equivalent 
commutation relations to half of an ebit. We can therefore 
use one ebit per two qubits to encode this code. 

In a later section, we show how to devise encoding and 
decoding circuits, beginning from k information qubits, c ebits, 
and a ancilla qubits per frame. Before that, however, we 
present a procedure to reduce the commutation relations of 
any check matrix to those of ebits and ancilla qubits. 

B. The Polynomial Symplectic Gram-Schmidt Orthogonaliza- 
tion Procedure for General Codes 

We detail a polynomial version of the symplectic Gram- 
Schmidt orthogonalization procedure in this section. It is a 
generalization of the procedure we developed for the above 
example. Before detailing the algorithm, we first prove a 
lemma that shows how to determine the shifted symplectic 
product matrix for an Z-expanded version of the generators 
by starting from the shifted symplectic product matrix of the 
original generators. 



1) The Shifted Symplectic Product Matrix for an l- 
Expanded Code: Suppose the shifted symplectic product ma- 
trix Q (D) of a given check matrix H (D) is as follows: 

Q(D)=Z (D) X T (D- 1 ) + X (D) Z T (D' 1 ) . 

Lemma 1: The shifted symplectic product matrix 17/ (D) of 
the ^-expanded check matrix Hi (D) is as follows: 



where the second line uses the invariance of the flooring oper- 
ation with respect to matrix transposition and time reversal, 
the third line expands the matrix multiplications using the 
matrix C[ (D) defined above, and the last line uses the matrix 
generalization of the multiplication property defined in ( fTOb . 
Our final step is to use the invariance of the flooring operation 
with respect to matrix addition: 



n, (D) = [ Rl (d^ 1 ) Q (d^ 1 ) Rf (D-V l )\ , ^ ( D ) = X l ( D ) Z f ( D ~ l ) + Z l ( D ) X T ( D 



where the flooring operation |_-J nulls the coefficients of any 
fractional power of D. 

Proof: Consider that the "X" matrix Xi (D) of the l- 
expanded check matrix Hi (D) is as follows: 



Xi 



i (D) = [Ri (d 1 ' 1 ) X (d 1 ' 1 ) C'i (D 1 / 1 ) 



where 



C, 



I (d 1/1 ) 



D° D 



D-C-D ] , 



and each diagonal D matrix is n x n-dimensional. It is also 
then true that the "Z" matrix of Hi (D) is as follows: 

Zi (D) = U (d 1 ' 1 ) Z (D 1 ' 1 ) C[ (d 1 ' 1 ) 

The matrix transpose operation, the time reversal operation 
(substituting D^ 1 for D), and matrix addition are all invariant 
under the flooring operation |_-J for arbitrary matrices M (D) 
<m&N(D): 



M D 



,i/i 



M T [D 1 ' 1 
M (d~ 1/1 
f N (D 1/l 



M (D 1 ' 1 

M [D 1 ' 1 



D=D~ 1 



N (D 1 ' 1 



Additionally, the following property holds for two arbitrary 
binary polynomials / (D) and g (D): 



[f (D 1 / 1 ) g (D" 1 

= J2 \p-*/ l f (D 1 ' 1 ) I [n^g^D 1 /' 



(10) 



Now consider the product Xi (D) Zj (D^ 1 ) 
Xi(D)Z?(D-') 
= Ri (d 1 ' 1 ) X (D 1 ' 1 ) C[ (d 1 ' 1 ) 

Ri^l^Z^D^Cl^D 1 ' 1 )^ 

Ri {d 1 ' 1 ) X (D 1 / 1 ) C'i (d 1 / 

l-l 



D=D- 1 



C'i 1 [D- 1 ' 1 ) Z T [D-^^Rj (D- 1 ' 1 



\ 

= £ [d-^ 1 Ri (d 1 ' 1 ) X (d 1 ' 1 

i=0 

\D l l l Z T (p- 1 ' 1 ) Rj [D- 1 ' 1 ) 

= Ri (D 1 / 1 ) X (d 1 ' 1 ) Z t (d- 1 ' 1 ) Rf (d- 1 ' 1 



R, D^ l )X D l / l )Z T D- l ' l )R 



Ri (D 1 ' 1 ) Z (D 1 / 1 ) X t (D- 1 / 1 ) R\ 

Ri (D 1 / 1 ) X (D 1 / 1 ) Z T {D- 1 ' 1 
+Ri(D l / l )Z{D l l l )X T (D 



-i/l\ 



T 

RJ{ 

Rf 



D -i/i 

D-i/i 

D -i/i 



r, D^n (d 1 i i \rT (d- 1 ' 1 



2) The Gram-Schmidt Procedure: We now present the poly- 
nomial symplectic Gram-Schmidt orthogonalization procedure 
that reduces the commutation relations of a given set of 
convolutional generators to have the standard form in (f3]). 

Consider the following n — k x 2n-dimensional quantum 
check matrix H (D): 

H(D)=[ Z(D) | X(D) ]. 

Label each row as hi (D) = [ Zi (D) I Xi (D) ] for all i 6 
{l,...,n-k}. 

We state the Gram-Schmidt procedure in terms of its effect 
on the shifted symplectic product matrix. It is easier to see how 
the algorithm proceeds by observing the shifted symplectic 
product matrix rather than by tracking the generators in the 
check matrix. 

Let I denote the amount by which we expand the check 
matrix H (D). Suppose first that I = 1 (we do not expand 
the check matrix). Let us say that the check matrix has r 
generators (we take a number different from n — k because 
the number of generators may not be equal to n — k for future 
iterations). There are the following possibilities: 

1) There is a generator hi (D) such that (hi hj) (D) = 
for all j 6 {1, ...,r}. In this case, the generator is 
already decoupled from all the others and corresponds to 
an ancilla qubit because an ancilla and it share the same 
commutation relations. Swap hi (D) to be the first row 
of the matrix so that it is hi (D). The shifted symplectic 
product matrix then has the form: 

•• 






h 





h-2.r 



r,2 



= [o]e 



i2 : 2 



*2,r 



where [0] is the one-element zero matrix and we use the 
shorthand hij = (hi hj) (D). We remove generator 
hi (D) from check matrix H (D) and continue to step 
two below for the remaining generators in the matrix. 
2) There are two generators hi (D) and hj (D) for 
which (hi hj) (D) = D m for some integer m and 
(hi hi) (D) = (hj hj) (D) = 0. In this case these 



generators correspond exactly to an ebit. Multiply gen- 
erator hj (D) by D m . This row operation has the effect 
of delaying (or advancing) the generator by an amount 
m and changes the shifted symplectic product to be 
(hi hj) (D) = 1. These two generators considered by 
themselves now have the commutation relations of half 
of an ebit. Swap the generators hi (D) and hj (D) to be 
the first and second respective rows of the check matrix 
H (D). Call them hi (D) and h 2 (D) respectively. The 
shifted symplectic product matrix is then as follows: 






1 


hl,3 ■ 


hl, r 


1 





tl2,3 ■ 


tl2,r 


^3,1 


/l3,2 


h3,3 ■ 


^3,r 



Ily 



'r,3 



We use the following row operations to decouple the 
other generators from these two generators: 



h'i (D) = hi (D) + {hi h 2 ) (D) ■ h Y (D) 
+ (hi hi) (D) ■ h 2 (D) for all i e {3, 



"}• 



The shifted symplectic product matrix becomes as fol- 
lows under these row operations: 






1 


• 





1 





• 











^3,3 • 


• K 



h 



r.3 



J( 



'3,3 



l r,3 



• h 



3,r 



■ hi. 



The first two generators are now decoupled from the 
other generators and have the commutation relations of 
half of an ebit. We remove the first two generators from 
the check matrix so that it now consists of generators 
h' 3 (D), . . ., h' r (D). We check to see if the conditions 
at the beginning of the previous step or this step hold 
for any other generators. If so, repeat the previous step 
or this step on the remaining generators. If not, see if 
the conditions for step three hold. 

3) There are two generators hi (D) and hj (D) for 
which (fn hi) (D) = (hj hj) (D) = but 
(hi hj) (D) ■£ D m for some m and (hi hj) (D) ^ 
0. Multiply generator hj (D) by 1/ (hj hi) (D). 
Generator hj (D) becomes infinite weight because 
(hj hi) (D) is a polynomial with two or more powers 
of D with non-zero coefficients. Now the shifted sym- 
plectic product relations are as follows: (hi hi) (D) = 
(hj hj) (D) = and (hi hj) (D) = 1. We handle 
this case as we did the previous case after the two 
generators there had the commutation relations of half 
of an ebit. 

4) None of these conditions hold. In this case, we stop this 
iteration of the algorithm and expand the check matrix 
by the next factor I := I + 1 and repeat the above steps. 

We have not proven that this procedure converges on all 
codes; however, it does converge on all the codes we have 
tried. We conjecture that this procedure converges for all 



codes, but even if this is true, in principle it might require 
expansion to a large number of generators. A simple and 
practical convergence condition is as follows. In practice, 
convolutional codes do not act on an infinite stream of qubits 
but instead act on a finite number of qubits. It may be that there 
are codes for which this procedure either does not converge or 
must be expanded until the frame size of the expanded code 
exceeds the number of qubits that the code acts on. In this 
case, we would not employ this procedure, and instead would 
treat the code as a block code, where we could employ the 
methods from Refs. Il23l . (24| f° r encoding and decoding. It is 
unclear if this practical convergence condition will ever really 
be necessary. This procedure does converge for a large number 
of useful codes, so that the frame size of the expanded code 
is much less than the number of qubits that the code acts on, 
and we have not found an example where this procedure fails. 

VI. Encoding and Decoding 

Entanglement-Assisted Quantum Convolutional 

Codes 

This section proves the main theorem of this paper. The 
theorem assumes that we have already processed an arbitrary 
check matrix with the Gram-Schmidt algorithm and that the 
shifted symplectic product matrix corresponding to the pro- 
cessed check matrix has the standard form in (01. The theorem 
shows how to encode a set of information qubits, ancilla 
qubits, and halves of ebits into a code that has equivalent error- 
correcting properties to those of a desired set of convolutional 
generators. 

The theorem uses both finite-depth and infinite-depth op- 
erations in the encoding circuit and finite-depth operations 
in the decoding circuit. The finite-depth property of the 
operations in the decoding circuit guarantees that catastrophic 
error propagation does not occur when decoding the encoded 
qubit stream. We do not review finite-depth or infinite-depth 
operations here but instead refer the reader to our previous 
article ifTTl that has a detailed discussion of these operations. 

Theorem 1: Suppose we have a set of quantum convolu- 
tional generators in the n~ kx 2n-dimensional matrix H (D) 
where 

H(D)=[ Z(D) | X(D) ]. 

Its shifted symplectic product matrix il (D) is as follows: 

Cl(D) = Z (D) X T (D- 1 ) + X (D) Z T (D- 1 ) . 

Suppose check matrix H (D) is the matrix resulting from 
processing with the polynomial symplectic Gram-Schmidt or- 
thogonalization procedure. Therefore, Vt (D) has the standard 
form in (0 with parameters c and a = n — k — 2c. Then 
there exists an online encoding circuit for the code that uses 
finite-depth and infinite-depth operations in the shift-invariant 
Clifford group IfTTl and there exists an online decoding cir- 
cuit for the code that uses finite-depth operations. The code 
encodes k + c information qubits per frame with the help of 
c ebits and a = n — k — 2c ancilla qubits. 

Proof: We prove the theorem by giving an algorithm to 
compute both the encoding circuit and the decoding circuit. 
The shifted symplectic product matrix D, (D) for check matrix 



H (D) is in standard form so that the first 2 c rows have the 
commutation relations of c halves of ebits and the last a rows 
have the commutation relations of a ancilla qubits. Perform 
the algorithm outlined in Refs. ||4], on the last a generators 
that correspond to the ancilla qubits. The algorithm uses finite- 
depth CNOT gates, Hadamard gates, and phase gates. The 
resulting check matrix has the following form: 



Z" (D) 
T(D) 



Z'(D) 




X'(D) 




(11) 



where the matrix Z" (D) and the null matrix at the top left 
of the "X" matrix each have dimension 2c x a, the matrices 
Z' (D) and X' (D) each have dimension 2c x n — a, and all 
matrices in the second set of rows each have a rows. The 
matrix V (D) may have entries that are rational polynomials. 
If so, replace each of these entries with a "1" so that the 
resulting matrix has the following form: 



Z" (D) 
I 



Z'(D) 




X'(D) 




This replacement is equivalent to taking a subcode of the 
original that has equivalent error-correcting properties and 
rate [3|. We can also think of the replacement merely as row 
operations with rational polynomials. We then perform row 
operations from the last a rows to the first 2c rows to obtain 
the following check matrix: 



Z'(D) 




X'(D) 




The shifted symplectic product matrix still has the standard 
form in (0 because these last row operations do not change 
its entries. 

We now focus exclusively on the first 2c rows because the 
previous steps decoupled the last a rows from the first 2c rows. 
Consider the following submatrix: 

H'(D) = [ Z(D) | X(D) }, 

where we have reset variable labels so that Z (D) = Z' (D) 
and X (D) = X' (D). Perform row permutations on the above 
matrix so that the shifted symplectic product matrix for H ' (D) 
changes from 



to become 



©■* 



/ 
/ 



(12) 



where each identity and null matrix in the above matrix is 
c x c-dimensional. 

We can employ the algorithm from Ref. |4| on the first c 
generators because the first c rows of the resulting check ma- 
trix form a commuting set (we do not use the row operations 
in that algorithm). The algorithm employs finite-depth CNOT 
gates, Hadamard gates, and phase gates and reduces the check 
matrix to have the following form: 





U(D) 




Z 2 (D) 



L(D) 

X X {D) 





X 2 (D) 



where L (D) is a c x c lower triangular matrix and U (D) is a 
cxc upper triangular matrix. The i th diagonal entry uu (D) of 
U (D) is equal to 1/lu (-D -1 ) where la (D) is the i th diagonal 
entry of L (D). This relationship holds because of the shifted 
symplectic relations in ( fT2l and because gates in the shift- 
invariant Clifford group do not affect the shifted symplectic 
relations. 

We now employ several row operations whose net effect is 
to preserve the shifted symplectic relations in (fT2l — we can 
therefore include them as a part of the original polynomial 
symplectic Gram-Schmidt orthogonalization procedure. Mul- 
tiply row i of the above check matrix by 1/lu (D) and multiply 
row i + c by 1/uu (D) for all i e {1, . . . , c}. Then use row 
operations to cancel all the off-diagonal entries in both L (D) 
and U (D). The resulting check matrix has the following form: 






I 





I Z! 2 {D) 


X[(D) 


X' 2 {D) 



where the primed matrices result from all the row opera- 
tions. One can check that the shifted symplectic relations of 
the above matrix are equivalent to those in (fT2l . Perform 
Hadamard gates on the first c qubits. The check matrix 
becomes 



/ 

X{(D) Z' 2 {D) 





1 X 2 (D) 



(13) 



We show how to encode the above matrix starting from 
c ebits and k + c information qubits. The following matrix 
stabilizes a set of c ebits: 



I I 






1 I 



(14) 



where each identity matrix is c x c and the last column of zeros 
in each matrix is ex (k + c). The receiver Bob possesses the 
first c qubits and the sender Alice possesses the last k + 2c 
qubits. The following matrix is the "information-qubit" matrix 

nn: 



where each identity matrix is (k + c) x (k + c) and each 
column of zeros is (k + c) x c. It is important to track the 
information-qubit matrix throughout encoding and decoding so 
that we can determine at the end of the process if we have truly 
decoded the information qubits. Perform finite-depth CNOT 
operations from the first c ebits to the last k + c qubits to 
encode the numerators of the entries in matrix Z 2 (D). Let 
Z 2 N (D) denote this matrix of the numerators of the entries 
in Z 2 (D). The stabilizer matrix becomes 



I I 








Z' 2 ,n (D) 



and the information-qubit matrix becomes 



Z» iN (D) I 






0/ 



where Z 2N (D) is the matrix that results on the "Z" side after 
performing the CNOT operations corresponding to the entries 



10 



in Z' 2 N (D). Perform Hadamard gates on the last k + c qubits. 
The stabilizer matrix becomes 



I 






ZLn (d) 



and the information-qubit matrix becomes 



Z'In ( d ) 




Let X 2 N (D) denote the matrix whose entries are the numer- 
ators of the entries in X 2 (D). Perform CNOT gates from the 
first c qubits to the last k+c qubits corresponding to the entries 
in X 2 N (D). The stabilizer matrix becomes 



/ / 

A (D) Z' 2N (D) 











I I X'(D) 



where A(D) = Z' 2N {D)X 2 ' N {D) and X 2 \ N (£>) is the 
matrix that results on the "Z" side after performing the CNOT 
operations on the "X" side corresponding to the entries in 
X 2 N (D). The information-qubit matrix becomes 






Z'l N {D) 


0/ 





X' 2 \ N {D) I 






Let r {D) be a diagonal matrix whose i th diagonal entry is the 
denominator of the i th row of Z 2 (D) and X 2 (D). We perform 
infinite-depth operations iPTD corresponding to the entries in 
r (D). The stabilizer matrix becomes 



/ 




r-^zr 1 ) o 

A(D) Z' 2tN {D) 




i r 




(D) 




Y' 

^2,N 


(°) . 


the information-qubit matrix becomes 






■ Z'l N {D)T-^(D- 
Xg iN (D)T- l (D- 








I " 







The above stabilizer matrix is equivalent to the desired one 
in (fT3] l by several row operations. We premultiply the first set 
of rows by V (D^ 1 ) and multiply the second set of rows by 
r _1 (D). We can also use the resulting identity matrix in the 
first set of rows to perform row operations from the first set of 
rows to the second set of rows to realize the matrix X[ (D). 
The operators that Bob would really measure need to have 
finite weight so he would measure the operators corresponding 
to the entries in the following stabilizer matrix: 



r 







I 
A(D) 








T(D) 





XLn (d) 



(15) 

We are done with the encoding algorithm. Alice begins with a 
set of ebits and performs the encoding operations detailed in 
(1141115b and then performs the finite-depth operations detailed 
in (II lti 1 3b in reverse order. 

We now detail the steps of the decoding algorithm. Perform 
CNOT gates corresponding to the entries in X 2 N (D) from 
the first set of c qubits to the last set of k + c qubits. The 
stabilizer matrix becomes 



I 
B(D) 



A(D) 




ZLn (D) 





1 T(D) 



where B (D) = 
matrix becomes 



Z 2 N (D) X 2 ' N (D). The information-qubit 



Y" 

^2,N 



ZL tN (D)T-i(D-i) 

(D) X"(D)T-i(D-i) 








Perform Hadamard gates on the last set of k + c qubits. The 
stabilizer matrix becomes 



/ 
B(D) 



r- 1 (d- 1 ) 

A(D) 




T(D) 



and the information-qubit matrix becomes 







7" 



(D)r-i(D-i) 

X" (D)T^(D^) 




ZLn (D) 




0/ 



Perform CNOT gates from the first set of c qubits to the last 
set of k+c qubits. These CNOT gates correspond to the entries 
in Z 2 N (D). The stabilizer matrix becomes 



/ 

B(D) 



A(D) 





T(D) 



and the information-qubit matrix becomes 



ZLn(D) Zi N {D)T^(D^) I 
X"{D) X>'{D)T^(D^) 








Row operations from the first set of rows of the stabilizer to 
each of the two sets of rows in the information-qubit matrix 
reduce the information-qubit matrix to the following form: 



Then we perform the finite-depth operations detailed in (fTTT 
n~3T >. We have now finished the algorithm for the decoding 
circuit because the logical operators for the information qubits 
appear in their original form. ■ 



A. Discussion 

Similar practical issues arise in these circuits as we dis- 
cussed previously in Ref. ifTTI . Encoding circuits with infinite- 
depth operations operations are acceptable if we assume that 
noiseless encoding is possible. Otherwise, infinite-depth oper- 
ations could lead to catastrophic propagation of uncorrected 
errors. Noiseless encoding is difficult to achieve in practice 
but we may be able to come close to it by concatenation of 
codes at the encoder. 

There is a dichotomy of these codes similar to that in Ref. 
[ 11 1. Some of the codes may have a simpler form in which the 
encoding circuit requires finite-depth operations only. Ref. lfT2ll 
gives an example of a code with this structure. These codes 
fall within the first class of codes discussed in Ref. [11] and 
will be more useful in practice because they propagate errors 
in the encoding circuit to a finite number of qubits only. The 
remaining codes that do not have this structure fall within 
the second class of codes whose encoding circuits have both 
finite-depth and infinite-depth operations and whose decoding 
circuits have finite-depth operations only. 
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B. Importing Classical Convolutional Codes over GF (4) 

One benefit of the new entanglement-assisted quantum 
convolutional codes is that we can produce one from an 
arbitrary classical convolutional code over GF (4). The error- 
correcting properties of the classical convolutional code trans- 
late to the resulting quantum convolutional code. It is less 
clear how the rate translates because we use the expansion 
technique. We know that the term (2k — n) jn lower bounds 
the "entanglement-assisted" rate 1 1 1 1 where n and k are the 
parameters from the imported classical code. The rate should 
get a significant boost from entanglement — the rate boosts by 
the number of ebits that the code requires. 

The construction for importing an [n, k] classical convo- 
lutional code over GF (4) is as follows. Suppose the check 
matrix for the classical code is an n — k x n-dimensional 
matrix H (D) whose entries are polynomials over GF (4). 
We construct the quantum check matrix H (D) according to 
the following formula: 



H(D)= 1 



ujH(D) 

luH(D) 



where 7 denotes the isomorphism between elements of GF[A) 
and symplectic binary vectors that represent Pauli matrices. 
Specifically, 7~ 1 (ft) = u)h x + uih z where ft, is a symplectic 
binary vector and h x and h z denote its "X" and "Z" parts 
respectively. We use this construction for the example in the 
next section. 

VII. Example 

We take the convolutional generators from Ref. (lOJ as our 
example. Ref. [23] originally used these generators in a block 
code. We import the following classical convolutional code 
over GF(4): 



(•••|0000|lwl0|1101|0000|---) 



(16) 



We produce two quantum convolutional generators by mut- 
liplying the above generator by ui and ui and applying the 
following map: 



(17) 



where II denotes the group of Pauli matrices for one qubit. 
The resulting quantum convolutional generators are as follows: 



GF(4) 


n 





i 


U) 


X 


1 


Y 


W 


z 



{■■■\IIII\ 



>,\IIII\---) 



IIII\ZXZI\ZZIZ\IIII 
IIII\XYXI\XXIX\IIII 



{■■■\IIII\ 



\IIII\...). 



(18) 



These generators have the following representation in the 
polynomial formalism: 



l + D 





D 
1 



D 







l + D 



1 
l + D 




D 



The shifted symplectic product matrix Q (D) for the above 
code is as follows: 

fD- 1 D- 1 
D D + D- 1 ' 



n{D) = 



D 



The above matrix is not reducible to the standard form by any 
row operations. We therefore expand the code by a factor of 
two to give four generators with a frame size of eight. The "Z" 
matrix Z (D) of the two-expanded check matrix is as follows: 



Z(D) 



1 10 110 1 

10 

D D D 1 1 

10 



and the "X" matrix X (D) of the two-expanded check matrix 
is as follows: 



X(D) 



10 

1 110 110 1 
10 
DDODlllO 



The shifted symplectic product matrix il 2 (D) of the two- 
expanded check matrix is as follows: 



Sh(D) 



l + D- 1 D- 1 

1 l + D- 1 

l + D 1 

D l + D 



We proceed with the Gram-Schmidt procedure because this 
matrix satisfies its initial requirements. We swap generators 
two and three to be the first and second generators of the 
check matrix because they have the commutation relations of 
half of an ebit. The shifted symplectic product matrix becomes 






1 





l + D- 1 


1 





l + D 








l + D- 1 





D- 1 


+ D 





D 






Multiply generator two by 1 + D and add to generator four. 
Multiply generator one by l + D -1 and add to generator three. 
The shifted symplectic matrix becomes 



1 

1 

l + D- 



D 2 







D- 1 





D~ 



We finally divide generator four by l + D- 
symplectic product matrix then becomes 



D 1 and the shifted 



m 



so that it has the commutation relations of halves of two ebits. 
The check matrix resulting from these operations is as follows: 



H 2 (D) = [ Z 2 (D) J X 2 (D) ] 



(19) 
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where 



Zl (D) 



D 

1 D 

D 
1 

1 




Xl(D) = 



1 


1 
1 


1 

1 + D- 

i_ 

D 





1 

1 1 





1 1 



-D- 1 

-D- 1 


-D- 1 



D +D 
D 2 +D+l 

D +D 
D 2 +D+l 



D 2 +D 
D 2 +D+l 

D+l 
D 2 + D+l 

1 
D' 2 + D+l 

D+l 
D 2 + D+l 



D 

D 2 + D + l 

D 
D 2 + D + l 



D 

D 2 + D+l 

D 2 + D + l 

D 
D 2 + D+l 

1 
D 2 + D + l 





(We give the transposition of the above two matrices because 
they otherwise would not fit in the space above.) The error- 
correcting properties of the above check matrix are equivalent 
to the error-correcting properties of the original two genera- 
tors. 

This code sends six information qubits and consumes two 
ebits per eight channel uses. The rate pair for this code is 
therefore (3/4, 1/4). 

We can now apply the algorithm in Theorem[TJto determine 
the encoding and decoding circuits for this code. The encoding 
circuit begins from a set of two ebits and eight information 
qubits per frame with the following stabilizer matrix: 



H (D) = [ Z (D) | X (D) ] 



where 



Z (D) 



X {D) 



1010000000 
0101000000 
0000000000 
0000000000 

0000000000 
0000000000 
1010000000 
0101000000 



We label the eight qubits on the right side of each matrix 
above as 1, . . . , 8. We label Bob's two qubits on the left as 
Bl and B2. Perform the following finite-depth operations (in 
order from left to right and then top to bottom): 

C (1, 4, D + D 2 ) C (1, 5, 1 + D 2 ) C (1, 6, 1) C (1, 7, 1 + D) 
C (2,4, D)C (2,5,1 + D)C (2,6,1) 
H(i,...,8)C(l,A,D + D 2 + D 4 ) C (l, 5, D 2 ) 
C (1, 6, 1 + D) C (1, 7, D 2 ) C (2, 4, D + D 2 ) 

C (2, 5, 1 + D) C (2, 6, 1) C (2, 7,1 + D), 

where we use the notation C (qi, qi, f (D)) to represent a 
finite-depth CNOT gate from qubit one to qubit two that imple- 
ments the polynomial / (D), H (qi,. .. ,qj) is a sequence of 



Hadamard gates applied to qubits qi through Qj in each frame, 
P (q) is a phase gate applied to qubit q in each frame, and 
C (q, 1/f (D)) is an infinite-depth CNOT gate implementing 
the rational polynomial 1/f (D). Alice performs the following 
infinite-depth operations: 



H (1,2) CI, 



C 2, 



1 + D- 
1 



1 + L>- 1 + D- 



tD- 2 
H(l,2) 



She then finishes the encoding circuit with the following finite- 
depth operations: 

H (1, 2) C (2, 3, 1) C (2, 5, 1) C (2, 6, 1) C (2, 8, 1) (20) 

P(2)H(3,...,8)S(2,3) 

C(1,2,1)C(1,3,1)C(1,5,1)C(1,6,1)C(1,8,1)P(2). 

The code she encodes has equivalent error-correcting proper- 
ties to the check matrix in ( fT9] i. 

Bob performs the following operations in the decoding 
circuit. He first performs the operations in (|20| > in reverse order. 
He then performs the following finite-depth operations: 

C (Bl, 4,D + D 2 + D 4 ) C (Bl, 5, D 2 ) 

C (Bl, 6, 1 + D) C (Bl, 7, D 2 ) C (B2, 4, D + D 2 ) 

C (B2, 5, 1 + D) C (B2, 6, 1) C (B2, 7,1 + D) 

H(3,...,8) 

C (Bl, 4, D + D 2 ) C (Bl, 5, 1 + D 2 ) C (Bl, 6, 1) 

C (Bl,7,l + D)C (B2,A,D) 

C(B2,5,1 + D)C(B2,6,1). 

The information qubits then appear at the output of this online 
decoding circuit. 

VIII. Conclusion and Current Work 

There are several differences between the methods used for 
general, non-CSS codes discussed in this article and the CSS 
codes used in our previous article ifTTl . It was more straight- 
forward to determine how to use ebits efficiently in CSS 
entanglement-assisted quantum convolutional codes, but we 
have had to introduce the expansion technique in Section |IV] 
in order to determine how to use ebits efficiently for codes 
in this article. There was also no need for an explicit Gram- 
Schmidt orthogonalization procedure in Ref. IfTTl . The Smith 
algorithm implicitly produced the row operations necessary for 
symplectic orthogonalization. 

We do have some methods that do not require expansion of 
a check matrix or an explicit Gram-Schmidt procedure |30|, 
but these methods do not make efficient use of entanglement 
and have a lower rate of noiseless qubit channel simulation 
and higher rate of entanglement consumption than the codes 
discussed in this article. Nonetheless, we have determined 
ways to make these other codes more useful by encoding clas- 
sical information in the extra entanglement with a superdense- 
coding-like effect [31]. These other codes are grandfather 
codes in the sense of Refs. |fl~2), 1132) because they consume 
entanglement to send both quantum and classical information. 
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One negative implication of the expansion of a code is 
that the expanded code requires more qubits per frame. The 
expanded code then requires a larger buffer at both the 
sender's and receiver local stations. The increased buffer will 
be a concern right now because it is difficult to build large 
quantum memories. This issue will become less of a concern 
as quantum technology advances. The entanglement-inefficient 
codes mentioned in the previous paragraph have the advantage 
that they do not require expansion and thus require smaller 
buffers. It therefore should be of interest to find solutions 
in between the entanglement-efficient codes discussed in this 
article and the entanglement-inefficient codes discussed in the 
previous paragraph. 

Some outstanding issues remain. We have not proven the 
convergence of the polynomial symplectic Gram-Schmidt or- 
thogonalization procedure and have instead provided a prac- 
tical stopping condition. We have a conjecture for how to 
proceed with proving convergence. Suppose that we would 
like to construct a code consisting of one generator that does 
not commute with shifts of itself. We have found for many 
examples that the correct expansion factor for the generator is 
equal to the period of the inverse polynomial of the generator's 
shifted symplectic product. We do not have a proof that this 
factor is the correct one and we are not sure what the expansion 
factor should be when we would like to construct a code 
starting from more than one generator. The conjecture about 
optimal entanglement use in general entanglement-assisted 
quantum convolutional codes from Ref. |29l also remains an 
open question. The conjecture is that the optimal number of 
ebits required per frame is equal to the following expression: 

rank (Z (D) X T (L*" 1 ) + X (D) Z T (D' 1 )) /2, 

where Z (D) and X (D) are the respective "Z" and "X" 
matrices for a given set of quantum convolutional generators. 
It is clear that this formula holds after we have expanded 
the original set of generators. The proof technique follows 
from the proof technique outlined in Ref. [29|. But we are not 
sure how to apply this formula to an initial set of unexpanded 
generators. 

The techniques developed in this article represent a useful 
way for encoding quantum information. This framework will 
most likely lead to codes with reasonable performance but 
they will most likely not come close to achieving the quantum 
capacity of an entanglement-assisted channel. The next step 
should be to combine the theory in this article with Poulin et 
al.'s recent theory of quantum serial-turbo coding [20|. This 
combination might lead to entanglement-assisted quantum 
turbo codes that come close to achieving capacity. 

MMW acknowledges support from NSF Grant CCF- 
0545845, and TAB acknowledges support from NSF Grant 
CCF-0448658. 
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